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where  V  =  voltage,  n  =  constant  greater  than  unity,  W  =  thickness,  A  =  constant  and  Ea  = 
oxygen  ion  mobility  activation  energy.  The  V  and  T  dependences  of  this  expression  agree 
with  our  experimental  results  for  Z5U  devices;  the  V  dependence  is  not  followed  by  X7R. 
This  expression  agrees  with  the  empirical  results  reported  by  Minford,  and  by  Procopowicz 
and  Vaskas. 
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SUMMARY  OF  RESULTS 


The  studies  reported  here  are  divided  into  two  sections:  Leakage 

Current,  and  Intrinsic  Degradation.  The  former  includes  transient  behavior, 
current -voltage -temperature  behavior,  studies  of  charge  carriers,  and  their 
transport.  The  latter  includes  time  dependence  of  degradation,  other 
degradation  characteristics,  and  degradation  models.  These  items  are 
summarized  below  in  that  sequence. 

1.  Leakage  Current:  Transient  Behavior 

Below  about  130°C,  especially  at  low  voltages  (<10V),  polarization  currents 
dominate.  Under  DC  bias,  capacitance  values  (dielectric  constant)  for 
commercial  Z5U,  X7R  and  COG  devices  (and  all  other  ceramic  types  measured) 
decrease  as  the  logarithm  of  time.  This  results  in  polarization  current  with  a 
time  dependence  of  t~m  ,  where  m  is  nearly  unity  near  room  temperature, 
decreasing  to  near  zero  at  ~120°C,  and  not  a  strong  function  of  voltage.  This 
current  dominates  over  short  times,  below  120°C,  and  is  generally  not 
correlated  with  conduction  (transport)  current.  Due  to  hysteresis  effects,  the 
transport  current  for  these  materials  cannot  be  reliably  determined  as  the 
difference  betwen  polarization  (charging)  and  depolarization  (discharge) 
currents.  The  displacement  current  consists  of  extra  terms  due  to  time  and 
voltage  dependences  of  the  dielectric  constant. 

2.  Leakage  Current:  I-V-T  Behavior 

For  all  samples  measured,  ohmic  behavior  at  low  fields  (~10^V/cm)  is 
followed  by  superohmic,  the  most  common  type  being  a  power-law  (Vn)  with 
l<n<3.  One  possibility  for  the  V^/2  relation  commonly  seen  is  electron  emission 
from  electrode  asperities,  which  were  identified  in  cross  sections,  and  are 
responsible  for  such  currents  In  other  solids.  The  near-square  law  behavior 
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often  measured  for  degraded  devices  is  attributed  to  a  "dulling"  of  these 
electron  emission  points  by  localized  ion  movement. 

Activation  energies  decrease  with  voltage  for  all  types  measured  other  than 
X7R,  accompanied  by  super-ohmic  voltage  dependence.  This  is  indicative  of 
GB-controlled  transport,  or  of  a  voltage -dependent  hopping  potential  barrier. 

3.  Leakage  Current:  Type  of  Carrier 

It  was  confirmed  that  the  charge  carrier  for  current,  for  the  type  of  X7R 
devices  and  blanks  measured,  is  the  electron.  This  assertion  is  based  on  the 
following:  a)  Leakage  currents  under  accelerated  conditions  are  constant  with 
time  over  extended  periods  (several  months)  prior  to  degradation;  b) 
Thermoelectric  voltages  are  negative  and  stable;  c)  No  Galvanic  voltages  could 
be  detected  for  X7R  blanks;  theoretically  predicted  Galvanic  voltages  were 
measured  for  Y203:Zr02  pieces  in  the  same  apparatus;  d)  Space  charge  limited 
currents  are  attributed  to  electrons  injected  from  the  cathode. 

4.  Leakage  Current:  Carrier  Transport 

Hall  voltages  could  not  be  measured  for  X7R  pieces,  even  at  elevated 
temperatures  (~200°C),  due  to  low  carrier  mobility.  Thermoelectric 
measurements  confirmed  the  low  mobility  (<10_®cm2/Vsec)  ,  and  that  thermal 
activation  energy  is  dominated  by  mobility  for  both  new  and  reduced  pieces. 
Carrier  concentrations  are  essentially  temperature  independent.  These  facts 
are  consistent  with  small  polaron  hopping  transport. 

The  role  of  grain  boundary  impedance  was  sought  using  complex  impedance 
measurements.  Evidence  of  grain  and  GB  contributions  were  seen.  Grain 
boundary  modelling  indicates  the  importance  of  grain  curvature,  and  size 
homogeneity,  in  maintaining  high  ceramic  impedance,  if  GB  controlled. 
However,  reasons  are  given  why  a  strong  case  cannot  yet  be  made  for  GB 
control  of  transport  for  our  X7R  pieces. 
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5.  Intrinsic  Degradation:  Time  Dependence 

Most  degradation  measurements  were  made  on  commercial  Z5U  and  X7R 
devices.  The  most  common  form  of  current  increase  with  time  for  both  types, 
prior  to  failure,  is  exponential  (I~e/^) .  Power  law  behavior  was  also  seen  for 
X7R  at  shorter  times.  For  Z5U,  degradation  rate  constant  increases  as  a 
power  law  with  voltage  and  exponentially  with  temperature.  For  X7R,  voltage 
and  temperature  dependence  are  both  best  fit  by  exponentials. 

6.  Intrinsic  Degradation:  Other  Characteristics 

The  exponentially  increasing  current  (with  time)  is  accompanied  by  a 
linearly  decreasing  activation  energy.  Thermal  activation  energies  decrease 
(from  ~1.3eV)  with  degradation,  being  zero  for  severe  degradation.  This  is 
attributed  almost  totally  to  transport  (reduced  hopping  or  grain  boundary 
potential.)  Ceramic  conductivity  varies  exponentially  with  activation  energy, 
except  for  severely  degraded  cases.  The  current-voltage  power  law  exponent 
increases  with  degradation;  we  attribute  this  to  the  "smoothing"  of  the  cathode 
("virtual  cathode"  effect.) 

For  simulated  degradation  (reduction  in  hydrogen),  carrier  concentration 
and  mobility  both  Increase  by  orders  of  magnitude. 

Concurrently  with  degradation,  color  gradients  developed  across  Z5U 
layers.  A  polarity  reversal  "healing"  effect  was  measured  for  Z5U,  X7R  and 
thick  film  capacitors  during  degradation.  These  effects  are  Indicative  of  ionic 
movement . 

7.  Intrinsic  Degradation:  Models 

Two  models  were  derived  for  the  degradation  rate  constant,  /?.  These  are 
both  based  on  oxygen  diffusion  as  the  cause  of  degradation,  the  effects  being 
Increased  carrier  concentration  and  reduced  transport  potential  energy  barrier. 
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The  rate  constant  is  found  to  have  the  following  voltage  (V),  temperature  (T), 
and  thickness  (W)  dependences: 


P=  AV^e-Ea/kT 
W2 

where  n  is  greater  than  unity,  A=constant,  and  Ea  is  the  oxygen  ion  mobility 
activation  energy.  If  lifetime  is  defined  as  ,  the  same  equation  obtained 
empirically  by  Procopowicz  and  Vaskas,  and  by  Minford,  is  obtained. 
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1.  INTRODUCTION 


This  report  discusses  work  performed  on  ONR  contract  N0014-83-K-0168, 
"Intrinsic  Mechanisms  of  Multilayer  Ceramic  Capacitor  Failure"  over  the  period 
of  the  contract,  3/15/83  to  5/14/88.  The  overall  objectives  of  the  program  were 
to  study  fundamental  (intrinsic)  leakage  current  and  degradation  behavior  in 
MLC  capacitors  and  related  material,  and  to  model  such  behavior. 

A  variety  of  measurements  were  made  to  address  these  objectives.  Since 
degradation  and  failure  in  an  MLC  device  almost  always  means  increased  leakage 
current  leading  to  shunt  or  short-circuit  conditions,  these  measurements  were 
aimed  at  understanding  the  nature  of  leakage  current,  and  the  causes  for  its 
increase  with  time  under  voltage-temperature  stress.  Therefore,  the  technical 
material  in  this  report  is  divided  up  into  two  main  sections:  Leakage  Current 
(Section  2),  and  Intrinsic  Degradation  (Section  3.) 

The  presentation  below  does  not  go  in  any  particular  chronological  order, 

1. e.  the  order  in  which  the  work  was  performed.  An  attempt  is  made  to  present 
results  in  a  manner  to  give  the  clearest  picture  of  our  findings  related  to  MLC 
current  and  degradation. 

Also,  since  most  of  this  material  has  been  presented  elsewhere  in  greater 
detail  (reports,  publications,  etc.),  it  is  presented  here  in  a  more-or-less 
summary  form.  Please  refer  to  the  more  detailed  references,  where  noted. 

2.  LEAKAGE  CURRENT 

Leakage  current  is  defined  in  this  report  as  current  caused  by  transport 
of  true  charge  across  the  dielectric.  This  is  contrasted  to  displacement  current 
which  is  enhanced  by  polarization  phenomena,  for  the  DC  case.  For  low  fields 
the  total  current  density  J  can  be  expressed  as 


5 


J  =  ctE  +  z0dE  ldi  +  dPldt 


(1) 


where  <t,  E  and  P  are  conductivity,  electric  field  and  polarization.  The  first 
term  represents  leakage  current  and  the  others  displacement  current. 

Direct  voltages  were  used  for  most  of  our  measurements,  including 
accelerated  stress  on  MLC  devices .  AC  voltages  were  used  mainly  for  the 
impedance  measurements,  to  be  discussed  in  Section  3.  DC  voltages  were  used 
for  all  of  the  measurements  discussed  in  Section  2. 

2 . 1  Transient  Behavior 

For  DC  studies  of  degradation  processes,  we  are  interested  mainly  in 
leakage  (transport)  current.  For  a  high  resistance  dielectric  near  room 
temperature,  this  is  difficult  to  measure  over  short  to  moderate  time  periods 
since  it  is  masked  by  polarization  current.  Such  transient  behavior  was  of 
interest  to  us  for  several  reasons:  i)  We  had  to  be  sure  that  polarization 
currents  had  decayed  to  zero  in  order  to  determine  transport  current;  ii)  It 
could  be  of  significance  with  regard  to  low  frequency  impedance  measurements; 
ill)  Such  transients  for  different  materials  may  be  of  interest  unto  themselves. 

Since  polarization  current  results  from  dipole  reorientation  which  changes 
polarization  P,  this  will  also  be  evident  in  capacitance,  as  shown  in  Figure  1, 
for  a  commercial  BaTiOj-based  X7R.  On  semi-log  plots,  capacitance  often  varies 
linearly  with  the  logarithm  of  time,  i.e. 

C  =  C0-C1lnt  (2) 


where  C0  and  Cj  are  constants. 

Polarization  current  Ip  under  DC  bias  can  be  expressed  as 
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X7R ,  50V 


▲  T  =  50  °C 
■  T  = 100  °C 
•  T  =  125  °C 


C(nF ) 


Figure 


TIME  (sec) 


Decrease  in  capacitance  with  time  for  X7R  device  at  three 
temperatures;  50V  bias.  The  solid  lines  are  least-square  fits  to 
the  equation  C  =  Cq  -  Cj  In  t. 


Ip  =  dQ/dt  =  VdC/dt  =  VCj/t 


(3) 


This  relation  is  fairly  closely  followed  at  low  temperatures,  as  seen  In 
Figure  2.  Over  the  time  scales  of  these  figures  (about  1  hour)  transport 
current  is  not  evident. 

Figure  3  illustrates  the  onset  of  transport  current,  after  approximately 
one  hour,  even  at  elevated  temperature  (100°C,  20V,  X7R) .  In  principle, 
transport  current  is  the  difference  between  the  magnitudes  of  polarization  and 
depolarization  (i.e.  short  circuit)  currents.  However,  for  ferroelectrics , 
polarization  currents  are  so  large  and  hysteretic  as  to  make  this  approach 
unreliable  (Figure  4). 

At  elevated  temperatures  and  voltages,  current  for  BaTiOg- based  X7R 
devices  decreases  with  time  as 


I  =  I0t-m  (4) 

with  exponent  n  varying  weakly  with  voltage  (Figure  5) ,  and  more  strongly 
with  temperature  (Figure  6).  Near  room  temperature,  n  is  close  to  unity  for 
these  devices,  decreasing  roughly  linearly  with  temperature  at  higher 
temperatures . 

The  current-time  behavior  of  BaTiOg-based  Z5U  devices  is  similar  in  form 
to  that  of  X7R,  as  indicated  in  Figure  4.  Exponent  m  values  are  near  unity, 
decreasing  with  temperature,  and  relatively  independent  of  voltage. 

The  agreement  between  measured  polarization  current  and  that  estimated 
from  equation  3  is  Indicated  in  Table  1. 
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Figure  2.  X7R  leakage  current  after  application  of  voltage,  for  five  voltages. 

The  decrease  nearly  follows  an  inverse  time  dependence  at  room 
temperature . 
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Figure  4. 


This  figure  illustrates  two  points: 
Z5U  and  X7R  devices,  although 
Transport  current  cannot  be 
depolarization  currents . 


i)  Equation  4  is  satisfied  by 
with  different  exponents;  ii) 
deduced  from  polarization  - 
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Figure  5. 
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The  exponent  parameter  m  varies  weakly  with  voltage  (X7R  shown 
Z5U  similar) . 
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Table  1 


Comparison  of  measured  currents  and  those  deduced  from  the  relation 


!d  =  civ/t  (X7R  at  100V, 

Z5U  at  10 V) 

Type 

Ci(F) 

t(sec) 

CiV/t(A) 

ID(measured)  (A) 

X7R 

3xl0-10 

100 

3xl0-10 

3xl0-10 

500 

6xl0-n 

5.6xlO-n 

1000 

3X10-11 

2.8x10-11 

2000 

1.5x10“ H 

1.4x10-11 

Z5U 

2.84xl0-8 

100 

2.8xl0-9 

2.9xl0-9 

500 

5. 7x10- 10 

7x10-1° 

1000 

2 . 8x10-1° 

3x10-1° 

2000 

1.4xl0-10 

1.6x10-1° 

From  these 

and  similar 

results  it  is 

evident  that  for  high  resistance  (i.e 

un-degraded)  ceramic,  true  leakage  currents  cannot  be  reliably  measured  at 
temperatures  less  than  about  150°C  unless  one  waits  for  times  on  the  order  of 
an  hour  or  more,  especially  at  lower  voltages. 


2.2  I-V-T  Behavior 

Useful  information  concerning  charge  carrier  injection  and  transport,  and 
degradation,  can  be  discerned  from  dependence  of  leakage  current  on  voltage 
and  temperature.  This  has  been  discussed  in  detail  in  progress  reports  [1-3], 
journal  articles  [4-9],  and  references  contained  therein.  The  major  types  of 
DC  current  we  have  considered  are  summarized  in  Table  2  [1,4]. 

I-V  characteristics  were  measured  for  new  and  degraded  capacitors,  mainly 
the  X7R  type.  New  devices  were  measured  at  elevated  temperatures  (>130°C) 
so  as  to  obtain  stable  currents;  degraded  devices  could  be  measured  at  lower 
temperatures  due  to  their  larger  transport  currents. 


Ohmic 


V 


Space  charge - 

planar  electrodes  Vn,  n  =;  2 

Space  charge -point  V^/2 

electrodes 

Schottky  exp  JV 

Several  trends  were  apparent  from  these  measurements.  Ohmic  behavior 
dominates  at  low  voltage  (<1V).  At  higher  voltages,  higher  power  voltage 
behavior  invariably  ensues,  the  most  common  type  being  at  or  near  a  V^/2 
dependence.  For  degraded  devices,  Vn  behavior  with  n>2  was  evident,  with 
hysteresis  (discussed  in  more  detail  in  next  section.) 

We  presented  a  model  for  the  V^/2  current  based  on  electron  emission  into 
the  dielectric  from  hemispherical  tips  of  electrode  protuberances  [1J.  Another 
example  of  a  somewhat  similar  nature  that  has  been  reported  is  electron  emission 
from  rough  cathodes  in  MOS  capacitors  (10-12).  Cross  sectioning  of  MLC 
devices  indicates  a  large  number  of  potential  electron  emitters,  electrode 
asperities  where  the  electric  field  will  be  enhanced  due  to  metal  curvature,  and 
where  the  dielectric  thickness  is  reduced.  Examples  are  seen  in  Figure  7. 
The  enhanced  electric  field  that  will  exist  at  the  tip  of  such  an  emitter  must 
be  viewed  as  a  possible  source  of  degradation,  since  both  local  joule  heating 
(due  to  Increased  current  density)  and  oxygen  ion  and  vacancy  diffusion  (due 
to  the  larger  field  and  perhaps  temperature)  will  be  enhanced.  Such  ion 
movement  could  account  for  the  transition  from  3/2  to  quadratic  voltage 


dependence  that  occurs  with  degradation.  This  is  discussed  in  that  context 
in  Section  3. 

Thermal  activation  energies  (<I>)  were  also  of  interest.  For  new  devices 
and  similar  BaTiC>3  ceramic,  values  of  1.2-l,3eV  are  typical.  There  are  several 
mechanisms  to  which  such  an  activation  energy  for  transport  current  may  be 
attributed:  oxygen  ion  diffusion  [14],  polaron  hopping  transport  [15,16], 

grain  boundary  (GB)  barrier  transport  [17,18],  and  activation  of  a  deep  donor. 
For  reasons  delineated  in  the  following  section,  oxygen  ion  diffusion  and  deep 
donor  activation  are  ruled  out,  leaving  polaron  and/or  GB -dominated  transport. 

We  discovered  and  reported  two  very  striking  properties  of  the  thermal 
activation  energies  for  MLC  devices:  they  decrease  markedly  (eventually  to 
zero)  as  devices  degrade;  they  decrease  with  applied  voltage  for  barrier  layer 
and  NPO  types,  although  not  for  the  X7R  type  from  two  vendors.  The  decrease 
of  4)  with  voltage  is  accompanied  by  superohmic  current.  Examples  are  seen 
in  Figure  8. 

Similar  examples  of  voltage-dependent  decreases  in  4>,  with  accompanying 
increases  in  current,  are  well  known  for  polycrystalline  silicon  [19]  and  ZnO 
varistors  [20]  (Fig.  8a).  In  such  cases  this  is  attributed  totally  to  collapse 

of  the  GB  potential  barrier.  Why  this  should  occur  for  three  of  the  device 

types  shown  in  Fig.  8  and  not  for  the  fourth  (X7R)  is  not  clear.  We  have 
discussed  the  role  of  the  GB  in  several  publications  [4-6, 8, 9],  The  role  of 
the  GB  with  regard  to  electron  transport  is  reviewed  in  Section  2.4. 

2.3  Type  of  Carrier 

MLC  capacitor  degradation  is  characterized  by  an  increase  in  leakage 

current.  Two  key  questions  addressed  in  this  contract  were:  what  is  the 

leakage  current,  and  what  causes  it  to  increase  with  time?  Over  the  course 
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Figure  8.  Current  (I)  and  activation  energy  (EA)  dependence  on  voltage  for 
a)  varistor,  b)  barrier  layer,  c)  COG  and  d)  X7R  devices 


of  this  program,  we  have  measured  mainly  commercial  capacitors,  or  similar 
"blanks”,  mostly  BaTlOg-based  X7R.  The  results  discussed  below  pertain  to 
such  samples.  The  nature  of  the  charge  carrier  is  discussed  in  this  section; 
its  increase  with  time  (degradation)  in  Section  3. 

Low-field  electrical  conductivity  is  proportional  to  carrier  concentration 
and  mobility.  Possible  carriers  are  ions,  electrons  and  holes.  In 
BaTiOg-based  material,  conduction  has  been  attributed  to  all  of  these.  Certain 
carriers  may  predominate  over  given  temperature  ranges,  or  oxygen  partial 
pressures . 

We  have  verified  that  the  dominant  charge  carrier  in  the  type  of  X7R 
ceramic  we  have  measured  is  the  electron.  (The  fact  that  a  majority  of  these 
electrons  may  originate  from  the  movement  of  oxygen  is  discussed  in  Section 
3.)  This  assertion  is  based  on  several  pieces  of  evidence: 

1)  Leakage  currents  under  accelerated  conditions  are  constant  with  time, 
over  extensive  time  periods.  An  example  is  shown  in  Figure  9,  where  current 
is  seen  to  remain  flat  under  different  acceleration  conditions  over  a  time  of 
nearly  three  months.  If  oxygen  vacancies  and  conduction  electrons  are 
generated  according  to  the  relation 

0o-*02  +  V  +2e'  (5) 

then  it  is  not  expected  that  leakage  current  would  remain  constant  if  vacancies 
and  electrons  (both  conducting  species)  are  generated  by  oxygen  diffusion. 

2)  Negative  Seebeck  (thermoelectric)  coefficients  indicate  that  the  charge 
carrier  Is  negative  [4J.  Reproducibly  negative  and  stable  voltages  generated 
verify  that  the  charge  carrier  is  negative,  and  of  constant  concentration  over 
the  time  of  the  measurement. 
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400  V,  170°  C 


Figure  9.  Leakage  current  versus  time  for  an  X7R  under  accelerated  stress 
conditions,  prior  to  degradation. 


2 


Seebeck  measurements  were  also  used  to  measure  changes  In  carrier 
concentration  and  mobility  for  BaTiOg-based  X7R  blanks  reduced  in  hydrogen. 
Thermoelectric  coefficient  S  versus  temperature  curves  are  shown  in  Figure  10a, 
for  "new"  and  reduced  samples.  Carrier  concentrations  were  deduced  from  the 
relation  [21] 


S  = 


(6) 


where  k/q  =  8.G3  x  10“  mV/K,  and  N  and  n  are  hopping  site  and  charge 
carrier  (polaron)  concentrations  respectively.  Drift  mobilities  were  determined 
from  the  relation  <t  =  nq/i ,  n  being  measured  simultaneously.  Carrier 
concentration  and  mobility  values  are  shown  in  Figures  10b  and  10c.  These 
are  discussed  more  in  Section  2.4. 

3)  Galvanic  cell  measurements  on  X7R  blanks  (supplied  by  Corning) 
indicate  that,  for  this  material  over  the  temperature  and  oxygen  pressures 
used,  the  charge  carrier  is  not  the  oxygen  ion.  If  there  were  considerable 
movement  of  oxygen  and  oxygen  vacancies  in  BaTiOg -based  material  (so  as 
to  make  a  contribution  to  conductivity,  above  the  electronic  contribution),  this 
should  be  detectable  as  a  Galvanic  voltage . 

The  Galvanic  setup  and  experimental  results  are  shown  in  Figure  11,  and 
were  discussed  in  detail  elsewhere  [3,6],  We  attempted  to  measure  Galvanic 
voltages  for  six  different  BaTiOg-based  samples:  two  each  of  Bag  ggTiOg  and 
®al  01TIO3  (provided  by  MRL,  Penn  State),  and  two  of  the  Corning  X7R  blanks 
mentioned  earlier.  A  yttrla  modified  zlrconia  (9.4M°/»Y203)  disc  was  used  as  a 
known  oxygen-ion  conducting  reference. 

Experimental  results  are  shown  in  Figure  lib.  The  Y203:Zr02  sample 
generates  voltages  slightly  less  than  predicted  by  the  equation 
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Figure  10.  For  BaTiOg-based  material  witli  X7R 
specifications,  in  as-received  and 
degraded  states:  a)  Thermoelectric 
coefficient;  b)  carrier  concentration; 
c)  mobility,  all  as  functions  of 
temperature . 
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which  is  the  solid  line  of  Figure  lib.  These  voltages  may  be  slightly  low 
because  the  temperature  (700°C)  is  somewhat  low  for  Zr02  [22], 

All  of  the  BaTi03-based  samples  generated  essentially  zero  voltage  at  all 
temperatures  and  oxygen  pressures  used.  This  indicates  that  oxygen  diffusion 
is  not  significant  in  these  samples  over  500-700°C,  and  is  therefore  probably 
not  significant  at  lower  temperatures.  Thus  it  appears  that  the  exponential 
increase  in  conductivity  with  temperature  seen  for  such  samples  is  not  due  to 
enhanced  diffusion  of  oxygen  or  oxygen  vacancies  (even  though  such  diffusion 
may  be  responsible  for  degradation  or  a  conductivity  increase  with  time,  as 
discussed  below. ) 

4)  We  have  reported  [1,4],  and  there  have  been  many  other  reports 
[23-26],  of  space  charge  limited  currents  in  BaTiOg-  based  materials  and 
capacitors.  Thus,  even  though  ionic  SCLC  can  potentially  exist,  and  has  been 
reported  in  other  materials  [27],  we  do  not  feel  that  the  SCLC  reported  for 
BaT103  is  ionic  because  there  would  have  to  be  a  steady  supply  of  ions  supplied 
by  an  electrode  (cathode  for  0“)  for  a  constant  current  to  exist  over  many  days 
or  even  months.  The  SCLC  is  due  to  electrons  injected  from  the  cathode. 

The  above  four  points  seem  to  verify  that  the  charge  carrier  forming  the 
leakage  current,  for  the  BaT103-based  samples  we  have  measured,  is  the 
electron . 


2.4  Carrier  Transport 

We  attempted  early  in  the  program  to  directly  measure  carrier  mobility  n 
using  the  conventional  van  der  Pauw  to  technique  [  1  ] .  These  results  were 
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unsuccessful,  i.e.,  no  Hall  emf  could  be  measured  even  at  elevated 


temperatures.  The  mobility  sensitivity  of  our  setup  was  about  O.Olcm^/Vsec. 
Other  measurements  (Seebeck  and  estimates  from  SCLC  characteristics)  confirm 
that  ft  values  for  such  ceramic  were  much  less  than  O.Olcm^/Vsec. 

For  undegraded  (un-reduced)  samples,  mobilities  were  determined  to  be 
less  than  about  10_^cm2/Vsec.  This  is  consistent  with  a  small  polaron  hopping 
mechanism  for  electron  conduction  [16?,  although  the  values  are  much  less  than 
those  reported  by  Seuter  [28]. 

The  conductivity  a  can  be  expressed  as 


<r  =  qn/i 


>/'oG 


-d>n/kTe  -<ykT 


(8) 


where  4>n  and  <I>/t  are  carriei  conceni  ation  and  mobility  activation  energies. 

Thermal  activation  energies  deduced  from  figures  10  and  11  for  X7R  blanks 
are  shown  in  Table  3. 


Table  3 

Thermal  Activation  Energies  (eV) 


Carrier 

Sample 

Concentration 

Mobilitv 

Conductivitv 

New 

0.19 

1.03 

1.22 

Reduced 

~0 

0.34 

0.34 

The  conductivity  activation  energy  is  ~1.2eV  for  as  received  X7R  chips, 
which  is  close  to  that  of  a  new  X7R  capacitor  from  the  same  vendor  (Corning.) 
The  decrease  in  activation  energy  of  the  reduced  chip  (0.34eV)  also 
corresponds  to  the  similar  decrease  we  have  measured  for  degraded  capacitors 
(discussed  in  Sec.  3.) 
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The  temperature  dependence  of  electron  transport  In  these  materials 
resides  almost  entirely  in  the  mobility.  The  weaker  temperature  dependence 
of  carrier  concentration  indicates  near-complete  ionization  of  donors. 

This  is  consistent  with  polaron  hopping  transport  of  electrons.  The  total 
electron  concentration  is  essentially  constant;  the  electrons  must  overcome  a 
potential  barrier,  which  would  be  the  of  equation  8  [38]. 

The  possibility  of  transport  being  dominated  by  grain  boundaries  (GB) 
was  treated  in  some  detail,  for  two  major  reasons:  1)  It  is  know  that  GB  control 
transport  in  certain  other  poly  crystalline  materials,  namely  thin  films  [29], 
polycrystalline  silicon  [19],  barrier  layer  capacitors  [30],  ZnO  varistors  [20], 
and  BaTiOg  thermistors  [17].  2)  The  voltage  dependent  activation  energy  that 

we  reported  for  COG  capacitors  (although  not  for  X7R.) 

Several  types  of  GB  potentials  are  shown  in  Figure  12,  for  homogeneous 
and  Inhomogeneous  grains.  Currents  that  can  flow  across  a  GB  are  also 
illustrated  in  Figure  12.  For  Schottky  barrier  type  electron  emission  over  the 
GB  barrier  (with  no  voltage  drop  in  the  grain),  the  I-V  relation  is 


I  =  I0e  °B/kT(eqV/kT -1)  (9) 

where  <I»g  is  the  GB  barrier  height.  At  a  given  voltage  this  expression  can 
be  approximated  by 


I  =  !Ae'EA /KT 


GO) 


where  EA  is  an  activation  energy  approximately  equal  to  <J>g  If  conduction  is 
dominated  by  electron  mobility  n,  the  latter  can  be  expressed 

/t  =  /i0e_(I>B^T  (11) 
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a)  GB  potential  energy  diagrams,  homogeneous  grains,  (a)  Negative 
charge  trapped  at  GB.  APS  =  0;  Qce  =  <7<VC8  tb)  Normal  polarization 
discontinuity.  APV  =  Qp.  Nob  =  0.  (c)  Reduced  poly-St.  no  charge  at  or 
near  GB  (flat  bands).  AP v  =  . V- 8  =  0 


b) 


(b) 

GB  band  diagrams,  inhomogeneous  grains,  (a)  High-resistance  layer 


near  GB.  (b)  Low-resistancc  layer  (c)  Presence  of  second  material  or 


phase 


Figure  12.  Grain  boundary  potential  energy  diagrams  for  a)  homogeneous 
grains;  b)  inhomogeneous  grains,  c)  Illustration  of  current  types 
in  the  vaclnity  of  a  GB. 
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Highly  super-ohmic  I-V  behavior  is  accompanied  by  decreasing  activation 
energy  (Figures  8a,  8b,  8c  for  varistor,  GBBL  capacitor  and  COG  capacitor, 
respectively).  The  weaker  superohmlc  behavior  of  BaTiOg-based  X7R  devices 
(^V1-5)  is  not  accompanied  by  a  decrease  in  (Fig.  8d)  We  reviewed  the  theory 
of  GB  transport,  and  voltage  dependent  GB  impedance  [5].  One  complication 
for  high  resistance  capacitor  ceramic  is  the  high  grain  resistance  (unlike  the 
GBBL  material,  or  any  of  the  others  mentioned  earlier) .  Even  so,  if  a  fraction 
of  the  applied  voltage  is  dropped  across  the  GB,  its  potential  barrier  should 
eventually  collapse,  resulting  in  superohmlc  current. 

Therefore,  the  question  to  be  addressed  was:  How  is  GB  impedance 

measured?  One  approach  is  to  place  micro-probes  directly  across  a  GB.  This 
is  feasible  for  large  grained  material,  especially  if  grain  resistance  is  small. 
However  for  smaller  grained  X7R  material,  with  high  grain  resistance,  the 
problems  of  small  size  and  surface  conductance  arise.  A  second  approach, 
which  we  pursued,  is  that  of  complex  impedance  measurements  of  the  type 
originally  reported  by  Baurle  [31].  The  reader  is  referred  to  the  references 
for  details  [3,9,31]. 

In  summary,  if  the  impedance  contributions  of  a  given  solid  can  be 
represented  by  lumped  R-C  sections,  these  can  be  discerned  on  a  complex 
impedance  plot  if  the  RC  time  constants  are  distinct.  The  complex  impedance 
Z  is  defined  as 


Z  =  Z'  +  j  Z"  (12) 

where  Z'  and  Z"  are  resistance  and  (capacitive)  reactance  respectively,  which 
are  both  functions  of  frequency  and  the  R-C  elements  representing  the  solid. 
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For  our  X7R  samples,  the  only  possible  RC  contributions  were  assumed 
to  be  from  the  GB  and  grains,  since  resistivity  was  not  a  monotonic  function 
of  sample  thickness,  thus  eliminating  contact  resistance  (see  Table  4) 

Table  4 

Resistivities  for  X7R  Chips  of  Different  Thickness 

Chip  Thickness  Resistivity 

(mm)  (12-cm) 


0.14 

1.59  x  106 

0.21 

2.31 

0.41 

1.77 

0.41 

1.92 

0.69 

1.90 

1.12 

2.17 

An  impedance  plot  (Z"  vs  Z')  is  shown  in  Figure  13a  for  a  commercial 
varistor,  across  the  frequency  range  5Hz  -  13MHz,  under  different  voltage 
biases.  Only  one  semicircle  is  evident  at  a  given  bias  because  only  the  GB 
contributes  significantly  to  impedance;  the  circle  radius  reduces  with  voltage 
because  of  GB  collapse. 

Many  impedance  plots  were  made  for  BaTiOg-based  samples;  all  showed 
evidence  of  more  than  one  impedance  contribution,  and  none  exhibited  a  single 
semicircular  structure.  An  example  is  shown  in  Figure  13b,  at  bias  voltages 
of  zero,  50V  and  100V.  Two  overlapping  semicircles  were  always  evident  for 
Corning  X7R  samples  of  this  type.  The  semicircles  are  offset  upward  from  the 
real  axis  (centers  indicated  by  crosses),  this  effect  attributed  to  local 
inhomogeneity  and  dispersion.  As  pointed  out  by  Kleitz  et  al.  [32],  the  offset 
Is  usually  small  for  intragrain  semicircles  (i.e.  the  grain  contribution.) 
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Figure  13.  Complex  Impedance  plots  for  a)  commercial  varistor,  and  b)  X7R 
ceramic,  at  different  voltage  biases.  The  key  for  the  X7R  Is-  #1 
zero  bias;  #2,  50V  bias;  #3,  100V  bias. 
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We  assigned  the  high  frequency  semicircles  (left  side)  to  the  grains,  and 
the  lower  frequency  ones  to  GB .  Resistance  values  obtained  from  this  data 
are  shown  in  Table  5,  for  different  biases,  for  three  samples.  It  is  seen  that 
the  GB  contribution  decreases  slightly  more  with  voltage,  which  was  generally 
the  case  for  all  samples  measured.  The  ratios  increase  with  voltage  bias  from 
two  percent  to  35  percent,  which  is  significant.  Such  changes  with  voltage 
are  consistent  with  superohmic  I-V  behavior  which  is  controlled  mostly  by  the 
GB. 

Table  5 

Rq/Rgb  Versus  Bias  Voltage 


- 

rg/rgb 

Bias  Voltage 
(V) 

1 

2 

3 

0 

0.58 

0.65 

0.46 

10 

0.73 

0.71 

0.51 

20 

0.70 

0.51 

35 

0.74 

0.69 

50 

0.72 

0.51 

100 

0.88 

0.61 

150 

0.47 

We 

also  reported  how 

grain 

size  could  be 

an  important  parameter 

with 

respect 

to  GB  impedance. 

GB  barrier  heights 

were  computed  versus 

grain 

diameter 

,  and  also  versus 

grain 

carrier  concentration  [9,33] .  Results  are 

shown  in  Figures  14  and  15.  Several  things  are  apparent  from  this  analysis, 
related  to  how  GB  can  Influence  ceramic  impedance: 

1)  Barrier  height  <J>  can  vary  over  orders  of  magnitude,  from  dominating 
(>leV)  to  non -important,  depending  on  grain  size  and  donor  density; 
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Grain  depletion 
region 


GRAIN  DIAMETER  (/tm) 


Figure  14.  Dependence  of  GB  barrier  height  on  grain  diameter  for  two  cases: 

ue  o  grain  depletion  only;  due  to  grain  depletion  plus  curvature 
Other  parameters  are  indicated  in  figure. 
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Figure  15.  Dependence  of  GB  barrier  height  on  donor  density  in  the  grain  for 
the  same  two  cases  as  in  Figure  14.  The  GB  barrier  height  peaks 
at  about  lO^cm-^ 


2)  Grain  curvature  can  cause  a  substantial  decrease  in  <I>  for  small  grains; 

3)  Resistance  homogeneity  throughout  the  material  could  depend  on  grain 
size  variation,  with  higher  curvature  grains  offering  reduced  resistance  paths; 

4)  Grain  size  variation  could  render  a  lumped  RC  model  for  the  GB 
Inappropriate;  a  distributed  model  would  be  neces-ary,  and  a  "distributed" 
impedance  plot  could  result.  This  could  also  account  for  some  of  the  impedance 
plot  offsets  which  we  measured. 

Thus,  even  though  these  results  seem  to  shed  some  light  on  the  role  of 
GB  in  high-K  ceramic,  their  exact  role  in  both  carrier  transport  and 
degradation  is  still  unclear  for  the  X7R  type,  the  main  one  we  measured.  There 
is  evidence  of  both  grain  and  GB  resistance  contributions  from  the  impedance 
plots.  However,  other  results  that  we  found  call  this  interpretation  into 
question;  namely  the  following: 

1)  Single  crystal  and  ceramic  BaTiOg  have  similar  activation  energies  E^, 
of  about  1.2-1. 3e V ; 

2)  The  two  apparent  semicircles  of  our  impedance  plots  have  similar 
activation  energies  of  about  1.3eV  [39]; 

3)  For  X7R  ceramic  and  devices,  values  are  essentially  independent 
of  voltage.  For  GB-  controlled  devices,  this  should  not  be  the  case. 

4)  BaTiC^-based  ceramic  resistance  has  been  reported  to  be  independent 
of  grain  size,  other  parameters  remaining  constant  [40], 

3  •  IMX&IMSiq  DEGRADATION 

Most  of  our  degradation  measurements  were  made  on  commercial  BaTiOg 
based  Z5U  and  X7R  capacitors,  under  accelerated  voltage  and  temperature 
conditions.  Other  pieces  of  similar  material,  without  internal  electrodes,  were 
also  studied.  Sections  3. 1-3.3  review  certain  electrical  characteristics  for 
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degraded  as  compared  to  "virgin"  samples.  Section  3.3  discusses  models  which 
can  account  for  these  characteristics. 

3.1  Time-Dependence 

By  far  the  most  common  type  of  time  dependence  seen  for  degradation  was 
exponential,  usually  following  some  initial  time  period  (typically  several 
minutes . )  Degradation  characteristics  for  X7R  devices  are  shown  in  Figs .  1G 
and  17.  After  an  initial  period,  current  varies  exponentially  with  ti’ne  until 
breakdown.  This  initial  transition  time  tQ  decreased  with  temperature  for 
measured  X7R  devices  as 


tQ  =  (4xl04  -  1 . 5x102T)  sec 


(13) 


between  about  200°C  and  280°C. 

Under  lower  stress,  leakage  currents  remain  constant  over  long  time 
periods  (Fig.  9).  We  have  also  reported  the  exponential  increase  for  thick  film 
capacitors  [34].  Z5U  devices  follow  a  similar  time  behavior  (Figure  18).  Thus, 
the  current-time  relation  for  capacitor  ceramic  can  be  represented  generlcally 
as  shown  in  Fig.  19,  where  the  time  scale  depends  on  voltage,  temperature, 
and  type  of  ceramic.  The  initial  decrease  in  current  was  discussed  earlier  (Sec. 
2.1). 

Such  capacitor  degradation  has  been  reported  previously  by  Prokopowicz 
and  Vaska  [35],  and  by  Minford  [36].  Relative  mean  lifetimes  were  described 
In  accordance  with  the  empirical  equation 


n 

T2 


=  (  )n  expC 
V1 


(14) 
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Figure  17.  Degradation  curves  for  X7R  devices  at  four  temperatures. 

Degradation  rate  parameter  p  was  determined  from  least  squares 
fits  to  such  curves.  For  these,  the  p  values  are  (xlO~^)  :/?j= 2.73; 
/?2=4.32;  p3  =6.91;  /?4= 34.7. 
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Figure  19.  Generic  current-time  curve  for  dielectrics. 
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where  tj  and  T2  are  device  lifetimes  at  voltages  Vj,V2  and  temperatures  Tj 
T2  and  Ea  is  an  activation  energy.  This  equation  has  shown  itself  useful  for 
numerous  MLC  capacitor  compositions,  at  least  at  lower  voltages  where  values 
of  n  tend  to  range  from  2.5  to  4  [37],  However,  this  was  an  empirical  equation. 
Over  the  course  of  our  project  we  have  developed  two  theoretical  derivations 
of  equation  14,  based  on  first  principles.  These  are  reviewed  in  Sec.  3.3. 

We  have  found  that  the  power  law  voltage  dependence  for  degradation  is 
followed  for  BaTK^-based  Z5U  and  thick  film  capacitors,  but  not  for  X7R,  which 
follow  an  exponential  behavior.  This  is  indicated  in  Figs.  20  and  21,  for  Z5U 
and  X7R  devices  respectively. 

The  exponential  temperature  dependence  of  degradation  rate  /?  was 
satisfied  by  Z5U  and  X7R  devices,  with  activation  energies  of  1.25  and  1.31eV 
respectively.  This  is  shown  in  Figures  22  and  23. 

3 . 2  Other  Degradation  Characteristics 

In  addition  to  the  direct  time  dependent  effects  noted  above,  several  other 
interesting  changes  in  the  properties  of  MLC  devices  occur  as  a  result  of 
accelerated  degradation.  These  effects  have  also  been  discussed  in  the 
publications,  which  should  be  consulted  for  more  details. 

a)  As  device  resistivity  decreases  during  degradation,  thermal  activation 
energy  decreases.  This  is  seen  in  Figure  24  for  an  X7R  device  aged  at  150°C, 
400V.  The  exponentially  increasing  current  is  consistent  with  a  linearly 
decreasing  activation  energy  since  I  ~  exp  (Ea/kT). 

Figure  24  pertains  to  a  single  device.  This  trend  applies  to  all  17  X7R 
devices  that  were  degraded  during  a  certain  test  early  in  the  program: 
resistivity  and  activation  decrease  monotonically  together.  This  is  shown  in 
Figure  25.  This  figure  is  consistent  with  Fig.  24  in  that,  for  the  upper  points, 


40 


DEGRADATION  RATE  (sec'1) 


VOLTAGE  (V) 


Figure  20.  Power  law  voltage  dependence  of  rate  parameter  ft  for  Z5U 
degradation,  at  three  temperature. 
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Figure  21.  Exponents!  voltage  dependence  of  rate  parameter  /?  for  X7R 
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Figure  22.  Arrhenius  plot  of  temperature  dependence  of  rate  constant  for  Z5U 
devices  at  different  voltages.  Activation  energy  =  1.25  +  O.leV. 
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Figure  23.  Arrhenius  plot  for  temperature  dependence  of  degradation  rate  for 
X7R  at  200  volts.  Activation  energy  in  region  2  is  1.31eV. 
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Figure  24.  Exponential  Increase  in  current  with  time  for  an  X7R  device  is 
accompanied  by  a  linearly  decreasing  activation  energy. 


5 


resistivity  increases  roughly  exponentially  with  activation  energy.  For  the 
lower  points  (resistivity  less  than  about  lO^flcm)  activation  energy  approaches 
zero. 

b)  As  ceramic  resistivity  decreases  and  leakage  current  increases,  the  I-V 
curve  becomes  more  super-ohmic.  For  the  device  of  Fig.  24,  the  voltage 
exponent  increased  from  1.25  to  1.65.  If  this  test  had  proceeded  longer,  the 
exponent  would  have  increased  somewhat  further,  and  then  dropped  to  unity 
as  the  device  approached  failure.  We  have  discussed  this  chronological 
evolution  of  I-V  characteristics  with  degradation  elsewhere  [4].  One  possible 
mechanism  Is  a  weakening  of  electron  emission  from  electrode  asperities  with 
time,  as  the  resistivity  of  the  ceramic  region  near  the  asperity  tip  decreases, 
due  to  accelerated  ionic  migration  In  the  enhanced  field  near  the  tip. 

c)  As  discussed  previously  for  "virgin"  samples,  the  activation  energy  for 
reduced  or  degraded  samples  resides  almost  entirely  in  the  mobility.  (Figure 
10,  and  Table  3).  It  is  seen  in  these  figures  that,  with  reduction  in  hydrogen, 
carrier  concentration  and  mobility  can  both  be  increased  by  orders  of 
magnitude.  According  to  the  models  discussed  later,  hydrogen  reduction 
represents  simulated  degradation  under  voltage-temperature  stress  because  both 
involve  loss  of  oxygen. 

d)  We  found  and  reported  color  gradients  that  developed  across  layers  of 
stressed  Z5U  devices  [1,4].  This  could  be  due  to  ionic  movement.  Some 
interesting  time -dependent  polarity -reversal  effects  were  also  seen  for  Z5U 
capacitors  (Fig.  26).  Exponential  or  near-exponential  current  increases  with 
time  are  evident.  In  Fig.  26a  it  is  seen  that  a  Z5U  device  whose  leakage  current 
Increases  by  two  orders  of  magnitude  heals  back  to  near  its  original  state  when 
the  voltage  is  removed  for  20  hours.  As  seen  in  Fig.  26b,  a  faster  healing 
can  be  produced  by  reversing  voltage  polarity.  Such  healing  effects  are 
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Figure  26.  Time  dependence  of  leakage  currents  for  Z5U  devices,  showing 
healing  effects,  a)  Resulting  f^-om  removal  of  bias,  b)  Resulting 
from  polarity  reversal.  * 
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indicative  of  ionic  movement  and  accumulation,  and  of  a  large  conductivity 
gradient  across  the  layer,  both  of  which  can  be  repaired  to  a  large  degree  by 
polarity-reversal  (with  similar  degradation  subsequently  occurring  in  the  other 
direction . ) 

It  is  of  interest  that  many  of  these  degradation  effects  may  be  so 
fundamental  as  to  apply  to  a  variety  of  dielectric  materials  (exponentially 
increasing  current,  near-linearly  decreasing  activation  energy, 
polarity-reversal  healing,  etc.)  As  an  example,  such  behavior  is  illustrated 
in  Fig.  27  for  BaTiOg-based,  K=500  thick  film  capacitors,  in  the  context  of 
another  program  [34],  Such  phenomena  can  be  accounted  for  by  ionic 
movement,  as  described  in  the  next  section. 

3.3  Models  for  Degradation 

We  have  developed  and  reported  two  models  to  account  for  the  degradation 
phenomena  that  have  been  seen.  The  first  one  is  based  on  oxygen  ion  and 
vacancy  migration,  leading  to  an  increase  in  conduction  electron  concentration, 
and  degradation  [3,6].  The  second  is  based  on  a  time-dependent  reduction 
of  grain  boundary  barrier  height,  resulting  in  increased  mobility,  and 
degradation  [34].  Since  the  role  of  grain  boundaries  in  high-K  ceramic  has 
not  been  conclusively  established,  and  since  the  grain  boundary  degradation 
model  also  depends  on  ionic  movement,  only  the  ionic  movement  model  is 
reviewed  here.  We  feel  that  such  a  model  might  pertain  directly  to  accelerated 
voltage-temperatures  degradation  phenomena  seen  in  high-K  ceramic,  and 
probably  to  many  other  dielectric  types  as  well. 

In  this  model,  there  are  assumed  to  be  two  sources  of  conduction 
electrons:  a  "background”  concentration  resulting  from  Impurities  and  defects 
in  the  as-made  material,  and  electrons  resulting  from  the  generation  of  oxygen 
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Activation  energies  versus  degradation  time  for  three  samples  aged 
under  indicated  conditions. 
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Figure  27.  For  K  =  500  thick  film  capacitors:  a)  Decrease  in  activation  energy 
with  degradation;  b)  Increase  in  leakage  current;  c)  Healing  by 
means  of  polarity  reversal. 
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vacancies  during  degradation.  The  generation  of  electrons  and  vacancies  is 
described  by  equation  (5),  with  the  electron  density  thus  becoming  equal  to 
2[V^J  (plus  the  background  concentration,  which  is  ultimately  neglected.) 
The  oxygen  vacancy  current  can  be  expressed  in  two  forms,  for  a  sample  of 
cross  section  A  and  thickness  W  [6].  These  are: 

Il0n  =  2qWA-^[V0-]  (15) 

and 

Iion  =  2qA,il^-[V--]  (1G) 

Equation  (15)  says  that  the  ionic  (0=)  current  is  proportional  to  the  time 
rate  of  increase  in  oxygen  vacancy  concentration.  This  assumes  that  none  of 
the  ions  that  leave  are  replenished  from  elsewhere. 

Equation  (16)  assumes  that  oxygen  ion  transport  in  the  solid  is  ohmic, 
where  V  is  voltage  and  ion  mobility.  This  assumption  gives  a  voltage 
power-law  dependence  for  degradation  that  is  too  low,  which  can  be  increased 
by  assuming  super- ionic  transport  for  the  ions,  which  results  in  a  voltage 
exponent  greater  than  unity. 

Equating  equations  (15)  and  (16)  we  get 


d[Vo  ] 

dt 


MjV 

W2 


[Vo] 


(17) 


It  is  now  assumed  in  the  degradation  regime  that  conduction  election 
concentration  nfe  =  2[V^  ].  Equation  (17)  is  then  integrated  to  yield 
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ne(t)  =  noe^t> 


(18) 


where  nQ  is  the  background  electron  concentration  and  n  Is  the  super-ohmic 
voltage  parameter  mentioned  earlier. 

The  empirical  result  reported  earlier  by  others  [35,36]  can  be  obtained 
by  defining  lifetime  r  as  equal  to  and  assuming  that  ionic  mobility  /q  is 

thermally  activated,  i.e.  /q  =  /i0  exp[  -Ea/kT] .  Lifetime  t  can  then  be  written 


T  =  /? 


-1 


w 


2  e  a 


En/kT 


n 


(19) 


The  ratio  of  lifetimes  for  devices  aged  under  different  voltage  and  temperature 
conditions  is 


(20) 


Equation  (20)  agrees  with  our  results  for  Z5U  devices  and  with  the 
empirical  lifetime  equations  reported  by  Procopowicz  and  Vaskas  [35],  and  by 
Minford  [36]. 
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